Cosmological perturbations in models of coupled dark energy 



o 



Sirichai ChongchitnarQ 
Oxford Astrophysics, Denys Wilkinson Building, Keble Road, Oxford, 0X1 3RH, United Kingdom. 

(Dated: October 2008) 

Models in which dark energy interacts with dark matter have been proposed in the literature to 
help explain why dark energy should only come to dominate in recent times. In this paper, we 
present a dynamical framework to calculate cosmological perturbations for a general quintessence 
potential and interaction term. Our formalism is built upon the powerful phase-space approach 
often used to analyse the dynamical attractors in the background. We obtain a set of coupled 
differential equations purely in terms of dimensionless, bounded variables and apply these equations 
to calculate perturbations in a number of scenarios. Interestingly, in the presence of dark-sector 
interactions, we find that dark energy perturbations do not redshift away at late times, but can 
cluster even on small scales. We also clarify the initial conditions for the perturbations in the dark 
' sector, showing that adiabaticity is no longer conserved in the presence of dark-sector interactions, 

I even on large scales. Some issues of instability in the perturbations are also discussed. 

^ . PACS numbers: 98.80.Cq 
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; I. INTRODUCTION 

Recent datasets from a range of cosmological experiments have convincingly suggested that roughly 25% of the 
I Universe is in the form of weakly interacting cold dark matter, whilst over 70% of the Universe is filled dark energy 
O with negative pressure that accelerates the cosmic expansion [3, [2) H, 3 • Despite these observational breakthroughs, 
-j-J I the theoretical origin of this dark sector remains poorly understood. Although most datasets can be explained by the 
^ ' ACDM model, in which dark energy is simply the cosmological constant, the data cannot yet rule out the possibility 
I— I , that dark energy may be dynamical, and that there may be undiscovered interactions in the dark sector (see [5] for a 
review). 

\ The simplest model of dark-sector interactions involves dark matter and dynamical dark energy in the form of 
I ' quintessence, a slowly-rolling scalar field evolving along a potential V{4>) [1, 0|- It is well-known that certain classes 
[ of V{4)) possess a desirable dynamical property such that field configurations from a wide range of initial conditions 
. evolve towards a common attractor. Furthermore, if there exists dark-sector interactions that allow dark matter to 
lO ■ continuously decay into dark energy, it may be possible to evoke such interactions to help explain the so-called "cosmic 
coincidence" problem, i.e. why dark energy only comes to dominate in recent times [sLlol. Il0|. 

The presence of dark-sector interactions would clearly affect the evolution of cosmological perturbations, which 
OO ! in turn affects the details of large-scale structure formation. Since specific cosmological details are clearly model- 
' dependent, it would be useful to have a general prescription with which cosmological perturbations could be calculated 
^ I given any model of quintessence and the dark-sector coupling. This is the primary aim of this paper. 

. There have been many investigations into the perturbations in coupled dark energy scenarios (see, for examples, 
/S\ [H [il [S El, [H, [H, El llll)- However, our approach is different as the calculation of perturbations will 
be entirely combined with the powerful phase-space formalism often used to analyse the dynamical attractors in the 
background [13, HT], [13, . As a result, the perturbation equations derived will be expressed in terms of dimensionless 
variables which are bounded, hence allowing a large range of initial conditions to be explored. We shall demonstrate 
the use of this technique with a few examples, which reveal interesting behaviour of the perturbations in the presence 
of dark-sector interactions. 

Some conventional views on cosmological perturbations are that matter density perturbations evolve adiabatically 
on large scales, and that dark-energy perturbations simply redshift away at late times. In the simplest models of 
quintessence, dark energy may cluster, but this is significant only on scales so large that it is irrelevant for astrophysical 
processes. We shall see that, in the presence of dark-sector interactions, all these statements are no longer true. 
Moreover, in some cases, instabilities can develop in the perturbations even though the background evolves smoothly. 
Throughout this work, we set c = 1 but keep the gravitational constant G explicit. 
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II. BACKGROUND EVOLUTION 



We assume a flat, homogeneous Friedmann-Robertson- Walker universe with metric 



(1) 



where r is the conformal time related to the cosmic time t hy dt = adr. We also assume that dark energy is a 
quintessence cj) with potential V{(j)), which is non-negative and at least twice differentiable in 0. 
The Friedmann equation is 



Pr + Pc + Pb + V + 



1 



(2) 



where k = VSttG = nip^ and pr, Pc, pt sue the energy densities of radiation, cold dark matter and baryons (an overdot 
indicates a derivative with respect to t unless stated otherwise) . The energy density and pressure of quintessence are 
given by 



1 



P4> 



1 -.^ 



V. 



The energy conservation equation 



Pi + 3Hpi{l + w.,) = 0, 



(3) 



(4) 



is assumed to hold for radiation (wr = 1/3) and baryons {wb — 0). The energy exchange in the dark sector can be 
represented by 



Pc + iHpc = Q, 
P0 + 3iJp0(l + w^) = -Q, 



(5) 
(6) 



where the coupling Q is a function of 0. In particular, if Q < 0, then dark matter instantaneously decays into dark 
energy. The reverse occurs when Q > 0. Nevertheless, it is important to note that the sign of Q does not by itself 
determine the overall effect of dark-matter-dark-energy conversion. For instance, it is possible for a model with Q < 
to end up with less dark-energy density at late times than a model with Q > 0, as we shall see later. 
Using Equations ([3]), Equation ([6]) becomes a modified Klein-Gordon equation 

4> + 3H:i, + ¥'{(/)) = -r, (7) 

where we define 

Q = r0. (8) 

The evolution of the various densities can be easily studied by introducing the dimensionless variables [13, HH, [l^l 



y 



K V 

i/V 3"' 



K Pr 



K Pc 



HV 3 '■ 



and, in addition, a dimensionless interaction variable 



The Friedmann equation then becomes the constraint 



kT 



(9) 



(10) 



x'^ + + z'^ + + v'^ = 1, 



(11) 



which means that the phase-space is simply the surface of the sphere S'^. 
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Next, differentiating each variable with respect to iV = Ina gives 



dN \2 dN 

where the A is the 'roll' parameter defined by 



d\nH 3 / , , 1 .\ 

dx /3, ^ { ^ d\nH\ , 

In - -^ + V2^^ -H' + ^j' ^''^ 

dy /3, fd\nH\ , , 

-Mo+-7^ ' (17) 



(18) 



Note that in this definition, = at the present time and decreases into the past. If A is constant, then V is just 
the exponential potential Voe^'^^'^. In general, any positive potential can be parametrized by a dynamical A(A^). 

In summary, given a potential V^(0), a coupling function Q{(\>) and initial conditions for {H,x,y, z,u,v), the set of 
coupled differential equations P^ - (fT71) completely specifies the evolution of the background cosmology. 



III. PERTURBATIONS IN A GENERAL MODEL 



While the dynamics of the background with dark-sector interactions have been well-investigated using phase- 
space techniques, the perturbations have so far been calculated by solving the perturbed field equations separately 
113, [13, HJj m, nil ■ This renders the two approaches somewhat disparate. 

In fact, it is possible to view quintessence models as trajectories in the phase-space {x,y,u,v,w), while simulta- 
neously studying the perturbations in the various energy densities along the trajectories. The main advantage of 
such an approach is that information about the attractors in the phase-space can be easily used to calculate the 
perturbations at (or near) various attractors. Initial conditions for the various densities are also easier to deal with 
since they correspond to finite regions in the phase-space. 

To demonstrate our technique, we consider a system with only dark matter and quintessence, and derive equations 
in the phase-space {x,y,u) governing the evolution of dark-matter density contrast, 5c, and that of the quintessence 
perturbation, 6(j), defined by 




Pc 



where an overbar indicates a background quantity. 



A. Perturbation equations 



As a starting point, we briefly follow the approach of Hwang and Noh [11] in their rigorous investigation of 'gauge- 
ready' forms of the perturbation equations for interacting fluids. 

Without fixing any gauge, the general scalar perturbations of the spatially flat FRW metric are given by 

ds^ = -a^{l + 2a)dT^ - 2a^(3,^dTdx' + ((1 + 2<y9)% + 2j^^j) dx'dx\ (20) 

where a, /9, (p, 7 are variables with small amplitudes. The components of the energy-momentum tensor are 

n = -p. Tl^-hp+p)v,,, t;:^p6;, (21) 
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where the total energy density p decomposes into the background value [pc + p^) and the perturbed part {5pc + Spcf,) 
and k is the wavenumber. The total pressure p and velocity v can be similarly decomposed. Overbars will now be 
omitted unless necessary. 

Hwang and Noh also derived perturbation equations in the 'comoving' gauge in which the dark matter peculiar 
velocity Vc vanishes. The variables 5c and 5(j) were shown to satisfy 



2HSc — ^PcSc 



-V't 



3H 



fc2 



)_d_ 

a2 dt 

r „ , 



r0 



Pc 

S<j) = (j)Sc + —5c - (1 + —)6T - -^S(j) 

Pc Pc Pc 

+ 2(0 + 3iJ0) — . 



dt \ Pc J 



(22) 



(23) 



Unless specified otherwise, the comoving gauge will be used for the rest of this paper. By changing the time variable 
to N and using the dimensionless variables ((9|)- p0)) and their evolution p^ - (fT7)) . the perturbation equations (|22p - ((23l) 
can be transformed into a dimensionless form 



dHc 



dm 



d^ (k(5 



dN'^ 



.dSc „ „dS^ „ dinSdi) , , 



(24) 
(25) 



which is a set of coupled second order differential equations, additionally sourced by the perturbation in the dark-sector 
interaction 57. The coefficients are given in terms of (x, u) and 7 by 



A 
B 
C 



dluH 



A = ^^^^^ + 2, 
dN ' 

2x 

E = 2V6a: + 4 



B ^ + ^ 



2a^^ + 4x7 + 4— rT^.T7 + VGXy J - 27^ 
dN dJy 



D = - 



2 rfln H 



din H ^ , 9 

4.x + 4 .„ X - 47 + vQXy 



dN 



87x2 



-7 H 

3 ' 3 dN 



1 d-f 



2 



F = 3y^X+4 



^^3dN-'^ 

^ dlnH dj , n 
57 + 47^ + -+47x^ 



2 27 



27AJ/ -3^/67 -2V67x--V67-^^^ ^ — 



27 Ay/e-f^x' 

7a ~^ 'i ^ 



dlnH 
dN 

iriy" + 
\/6x, 



^ + 3, 
aH 



37 + 27 



dlni? 



(i7 272X 
dN ^ u2 



2x 47^ 



2V6 



7A?/2 



Z? = 2V62:^ 



-V6(l + ^), 



87^x2 



(26) 



where d\nH/dN = —3(1 + x^ — 2/^)/2. We also define the curvature 



(27) 



These equations completely describe the perturbations along the trajectories in the phase-space {x,y,u). They are 
completely general and valid for any quintessence potential, parametrized by {A, rj}, and any interaction rate 7. Notice 
that most terms in the coefficients depend on 7, its derivatives, or its perturbation ^7. 
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B. Perturbed interaction 

It remains to show how 5^ may be calculated from a given interaction term Q which appears in the energy exchange 
equations dH])-®- 

First, we note that the introduction of Q amounts to a modification of the Einstein equation, and so one must be 
able to associate Q with some covariant entity, say, Q/^: which satisfies 

^.Tl'^'^" = Qf, (28) 
V,Tl^>^Q^^\ (29) 

Q'^^^^-Qi^^ (30) 

where the total energy-momentum tensor = T^*^^"^ + t'^^'^ satisfies the conservation ~ 0. The vector Q^^'' 

can be decomposed into the background and perturbed parts as [lH 

Q[,'^ - -a (q^'\i + a) + SQ^'^iy (31) 

Ol'^ = J^f- (32) 
Inserting these into ((28)) and using ([2T|) . we find 

Q^'^ = Q, a = — . (33) 

Pc 

The exact forms of SQ^'^'^ and J*-*^-* can be identified by perturbing the vector Qjf^ and comparing with Equations 
([3T|l - p3p . We shall demonstrate this procedure for a few examples in the next section. 

Next, recall that earlier we introduced F — Q/4>- I3y writing Q^''-' — — F(/)^^ and making the perturbations F T+6T 
and (j) ^ (j) + 6(1), Equations (l28])-(l33l) imply 

6r = (^Q(") - T{S4> - a</>)) /0, (34) 

a = (35) 

Pc 

Finally, having obtained F and ST, multiplying them by k/\/6H^ gives the dimensionless variables 7 and S'f, which 
appear in the main perturbation equations (|24|) - ((25)) . 

Note also that in the case without interaction (F = 0), we find a = 0, which is just the synchronous gauge condition, 
as observed in The synchronous and comoving gauges do not coincide if there is any dark-sector interaction. 

IV. APPLICATIONS 

We now apply Equations ([24 |) -([25 |) to calculate the perturbations in two models with simple interaction terms found 
in the literature. By setting 7 = 0, we also easily obtain perturbations in the case without interaction. 

In addition to the dark-matter density contrast 5c and the field perturbation d^i, we can also calculate the dark- 
energy density contrast = 5p^/p(f,, which quantifies the clustering of quintessence, as follows. The energy- 
momentum tensor of quintessence is 



(36) 



By perturbing the 00-componcnt and equating with r'"^^*^ = — (p^ + Spc/,), we find 



Sp^ = ^5^ - 4>'^a -f V'5(j,. (37) 
Dividing the above by p^ and using ([3]) and psp . we arrive at a useful expression for the dark-energy density contrast 



6d, = 



3(x2+2/2) 



V.2 



aN V 



(38) 



A. Initial conditions 



If there exists an attractor with a sufficiently large basin of attraction in the phase-space, then the evolution of the 
various background densities are generally insensitive to the initial conditions for (x, y, u). However, the evolution of 
the perturbations Sc and Scj) may not necessary be insensitive to the initial conditions in Sc and 6(j). These conditions 
determine the relative amplitudes of the various growing and decaying modes, and, therefore, must be carefully 
specified. 

In this paper, we shall impose the adiabatic initial condition, which specifies that the entropy perturbation between 
dark matter and quintessence and that intrinsic to the quintessence vanish at early times. Let S and T denote the 
relative and intrinsic entropy perturbations respectively. The adiabatic initial condition implies that 



Pc 



p<t> 



P4> P4> 

In terms of the dimensionless variables and field perturbations, the above conditions correspond to 



S = ^ 6c 

1 = ^ - 



_ - 272: 
V6(7 + 3a;)(K50) 



V6: 



d{KS(j)) 



dN 
d{KS(j)) 



3Xy K, 



dN 



2^67^ 



3A?/ K, 



(39) 
(40) 

(41) 
(42) 



Combining these conditions, we see that adiabaticity implies that the initial matter density contrast is proportional 
to the initial field perturbation, with 



^272; 



3u' 



k54>. 



(43) 



It is well-known that in the case with no interaction (7 = 0), adiabaticity is conserved on large-scales [23, l2a |. 
However, interactions in the dark sector can source isocurvature perturbations even for large-scale modes that are 
initially adiabatic. For example, using the adiabatic conditions, the expression for the derivative of iS is 



dS_ 

dt 



1 



k 



3%/6x \aH 



fil) 



(k(5(/)) 



3(7) - 



1 

3a; 



2x 
3^ 



(44) 



where 7(7) = 5(7) = whenever 7 = 0. Thus, we see that even in the large-scale limit, dS/dt ^ unless the 
interaction and its perturbation vanish. A similar expression for dX/dt can be found, albeit more complicated. 
Consequently, one cannot consistently insist that both S and dS/dt vanish initially, as was done in [isj . The key 
point here is that in the presence of dark-sector interactions, adiabaticity is generally no longer conserved on large 
scales. 



B. Example I: No interaction 



First, we check that the system (|24l) - ((25)) correctly reproduces well-known results in the case without interaction. 
Setting 7 = 57 = gives 



£5c 

dN^ 

d^ {KS(t)) 
dN"^ 



dlnH \ dSc 3 ,^ „ ^ dinSch) „ 9,, , 



dlnH 
dN 



d {K.5(j)) 
dN 



3r/y 



k 



(k(5(/)) = ^/&x 



dSc 
dN' 



(45) 
(46) 



Phase-space analyses [23, show that during the matter-dominated era, there exists a global attractor- the 
so-called scaling solution, whenever A > VS. During the scaling regime, we have 
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To see how the perturbations evolve during the scaling regime, we insert the coordinates of the scaling solution into 
P5)l - (|46p and seek solutions of the form 6c oc S(j) (x exp{pN){— qP). In the large-scale limit {k <C aH), the system can 
be solved analytically. The solutions are 



3 3 3 /24 
^=1' -2' -4±4Va^-^- 

These solutions agree with those found in (23 | . The first two modes are the well-known growing and decaying mode in 
the absence of quintessence. The quintessence induces two additional modes whose nature depends on A. Quintessence 
density contrast ([38]) also grows as at the scaling solution. 

Phase-space analyses also reveal another attractor for A < \/3, at which quintessence completely dominates the 
cosmic densities. In this regime, 



A / A2 d\nH A2 , ^ 

By proceeding as above, one finds that the resulting coupled system can also be solved analytically. In the large- 
scale limit, there is now one scale-invariant mode (p — 0), and three additional modes corresponding to the solutions 
of the cubic equation 



A^ 9A2 



4 



p+Y-3A^=0. (50) 



In the limit of a flat potential with A ^ (the so-called 'skater' models [30]) or potentials with very gentle slope, the 
perturbations consist of i) a constant (scale-invariant) mode Sc oa S(f> — constant, ii) a decaying mode Sc oc and iii) 
a decaying mode S(j) oc a~^. This behaviour of 5c during the dark-energy dominated epoch agrees with standard results 
(see e.g. [31|)- By setting A = 7 = in p8|) . we also find that the quintessence density contrast is exponentially 
damped towards zero in this limit. Thus, quintessence cannot cluster during a dark-energy dominated era. This, 
however, is not necessarily true if there are interactions in the dark-sector, as we shall see shortly. 

We now turn to the case with nonzero interaction. In each of the examples that follows, we summarise the 
background dynamics in terms of attractors in the phase-space. We then calculate the perturbations as the trajectories 
approach the various attractors. 



C. Example II: Q = f3Hpc 

The interaction of the form 

g = l3Hpc. (51) 

appeared in an earlier work of Lima et al. (3^ as a particle-creation mechanism during inflation, but has subsequently 
been studied in the context of dark energy by a number of authors ( [H, [H, [s^l , amongst others). It is not straight- 
forward to see what vector correctly gives rise to this form of Q. In a recent work, Valiviita et al. [l2] carried out 
an analysis using 



<^ (52) 



(c) 

where w/^ is the dark matter velocity 



u^^^ =a{-l-a,{vc + f3),i). (53) 



Physically, ([5^ means that the rate of interaction of dark matter is greatest along the direction of its velocity. 
By perturbing this vector using ([53]) and equating with ((3T|) . we deduce that 



5Q^^^ = Q5c- (54) 



In terms of the dimensionless variables, we have 
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Substituting all this into ([24)l - (l25p gives the coupled system 



dm 



2\/6a 



f d\nH \ ddc 3 2c 

/3(9 - 2/3) d In iJ _^ 5/3(3-/3) , Xy'^p{/3~3) u^f3^{l3~3) 



Pis - /3) 
\/6x 



d {K,5(j)) 
dN 



3y^X + 2\/6/3a 



/3 



dlnH 
dN 



dN 
+ 3 + /3 ( 1 

1 



2x2 



2x2 



2v^x3 



k 

X \ aH 



d (nScj)) 



dlnH 
dN 



3- 



dN 
3%/6A2/2 



3f]y- 



(- 



2x 



{n64>) = v^a;^ 



2x 



Note that by setting /3 = 0, we recover the system (|45I) - P5)) for the non-interacting case. 



At the scaling solution 



(57) 



(58) 



If A and /3 satisfy the conditions [2S 



<3, 



-2/3 < ( < 3 - /3, 



then the scaling solution exists, with 



3-^ 



y = 



6A2 ^ 3 



u = \ l- 



(3 - P? /3 
3A2 3 



dlnH _ f3-3 
dN ~ 2 



(59) 



(60) 



This is a global attractor and hence the growth rate of perturbations are insensitive to initial conditions in (x,y,u). 
When the phase-space is projected onto the x-y plane, the locus of the scaling solution is the hyperbola 



y - X - 



A 



_A2 
6 ■ 



At the scaling solution, the dark-energy density is given by 

= x^ + j/^ = 



2 , 2 (3-/3)2+/3A2 



3A2 



(61) 



(62) 



As mentioned earlier, it is not the sign of Q but, in this case, the expression ([62|) that determines the overall late-time 
effect of the coupling on the dark-energy density. Indeed, if Jl^j^^o is the dark-energy density in the case without 
coupling, then 



A2+/3- 




< f^</)|/3=0- 



(63) 



The effect of the interaction satisfying the first case in (I63p allows more dark energy to be produced at late times 
compared with the non-interacting case. This may be of some use in alleviating the coincidence problem. 

The effect of the coupling of the type (|5ip on the background is summarised in Figure [1] which shows the phase 
space projected onto the x-y plane. In the figure, we set A = \/5 and /3 ranging from —0.2 to 0.04 (right to left). 
These parameters satisfy the conditions ((59)) for the existence of the scaling solution. From a common set of initial 
conditions, trajectories are evolved towards the scaling solutions, which vary very slightly in position along the locus 
(pT|) . shown by the dashed line. Increasing the value of /3 effectively lengthens the path towards the attractor and 
pushes the trajectories to explore regions with smaller x. Physically, this means that the interaction with positive 
values of /3 effectively acts as an opposing force to the motion of the field. Similarly, the interaction with /3 < 
accelerates the field towards the scaling regime. 

To solve for the perturbations during scaling, one might try substituting the coordinates into ([57]) - ([55)1 . This, 
however, fails to yield a system of differential equations that can be solved exactly, due to the complicated dependence 
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FIG. 1: Trajectories in the x — y plane for the exponential potential V = Voe~"^'* with A — \/5 and interaction term of the form 
Q = PHpc. From left to right, the trajectories correspond to P = 0.04, 0.035, 0.02, 0, —0.1 and —0.2. The behaviour near the y 
axis seen in the left-most trajectory is discussed in the text. The dotted hyperbola is the locus of scaling solutions [Equation 
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FIG. 2: Evolution of large-scale perturbations over AA'" = 10 during the scaling regime for the exponential potential with 
A = \/5, and interaction term of the form Q = (3Hpc (with /3 = 0,±0.1). The panel on the left shows the evolution of the 
dark-matter density contrast 5c, while the panel on the right shows the quintessence density contrast \S,jy\, showing similar 
trends. Dashed lines indicate negative values of the perturbations. 



on the parameters /3 and A. It is nevertheless useful to integrate the equations numerically. In all numerical examples 
that follow, we impose the adiabatic conditions ([1T|) - (|^^ and set the initial amplitude Sc ~ 10^^ at the initial time 
Nq = —10. The remaining initial conditions can be obtained from the relations (|4T|) - (l43)l . 

Figure [2] shows the results of numerically integrating the dark-matter and quintessence density contrasts during the 
scaling regime of a large-scale mode (with k/aH > 10'^ in the range — 10 < < 0). The adiabatic condition favours 
a large initial amplitude of the decaying mode. The gradient of the negative slope is approximately —3/2 in all cases 
(see Equation (j48p ). The growing mode subsequently dominates. The gradient of the positive slope is approximately 
1 and increases slightly with negative values of f3. Similarly, positive values of /3 slows down the growth of the density 
contrasts. 

Dashed lines in the figure indicate negative values of the perturbations. Comparing the two panels in Figure [2l we 
see that at late times, clumping of dark energy is accompanied by dilution of dark matter (and vice versa). 



2. Instability near x = 



As we have seen, positive values of f3 work in opposition to the motion of the quintessence. When the quintessence 
slows down sufficiently so that x is close to zero, the background evolution must be handled with care, for \dx/dN\ oo 
(this stems from the standalone 7 term in Equation (|13p . with 7 = Pu^ /2x). This does not necessarily mean that the 
background is unstable, but it implies that any trajectory intersecting the y axis stays on the y axis. Putting x = 



10 




FIG. 3: Evolution of large-scale perturbations over AA'^ = 10 during the quintessence-dominated regime for the exponential 
potential with A = 0.01, and interaction term of the form Q = (3Hpc (with /3 = 0, ±0.1). The panel on the left shows the 
evolution of the dark-matter density contrast 5c, while the panel on the right shows the quintessence density contrast 
Dashed lines indicate negative values of the perturbations. Asymptotic forms of these curves are given in the text. Note 
especially that 5^ is no longer constant if /3 7^ 0. 



into (HH), y can be solved exactly as 

y=(Ce-3^' + l)-l/^ (64) 

with C constant. This means that whenever a trajectory intersects the y axis, it then travels vertically towards (0, 1), 
which is the quintessence-dominated attractor. This behaviour can be partially seen in the left-most trajectory in 
Figured] The numerical integration handles the trajectory so far as the intersection with the y axis, at which point 
it breaks down. 

Whilst the background remains well-defined when the field comes to a stop, the perturbations are not so well- 
behaved. Almost every coefficient in ([57|) -([55 1) blows up as a; ^ implying that neither 5c nor 5^ are stable at a 
linear level of perturbation. Hence, initial conditions in {x,y,u) must be tuned so that x never reaches 0. This sort 
of coupling is therefore not theoretically viable if quintessence is to alleviate the fine-tuning of initial conditions. This 
conclusion applies to any quintessence potential, and naturally extends to phcnomcnological interactions of the form 

M 

Q(xHipc + p^). (65) 

We note that a different type of instability in a similar model was identified in fT5| during the radiation-dominated 
era, in which the perturbations grow at extremely high rates. 



3. Flat limit 



For potentials with very gentle slope or flat 'skater' models, the perturbations for large-scale modes can usually be 
solved analytically. For these models, the quintessence-dominated attractor exists and is stable whenever 

< 3 - /?. (66) 

By solving ([571) - ([55)) and using the limit A ^ and k ^ aH, we find the usual scale-invariant mode 5c = constant, 
and 5c ~ e~^^ as before, plus two quintessence-induced modes. For /3 < 0, we find the decaying modes 

5c 5(1)^ e-'^^iCiCOseN + C2smeN), (67) 

where Ci's are constant and 9 = y/\5(3'^ + 30/3 + 9|. For (3 > 0, we find 

5c-- 5(1) e-'^^iCscoshON + C4 sinh 6'A^) . (68) 
The salient feature of these solutions is that, for /? > 0, the growth rate of the mode 

5c^5^^a-^+^ (69) 
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surpasses that of the scale- invariant mode (since the exponent is positive for /? > 0). Hence, Sc and 5^ could indeed 
grow during dark-energy-dominated era if /3 > 0. 

Figure [3] shows numerically evaluated perturbations for a few values of (3, illustrating features outlined above. For 
these models, we use potentials with A = 0.01, as A = would blow up the perturbation equations in the numerical 
program. The growing modes of the density contrasts when /? > are clearly seen in both panels. The evolution of 
the quintessence density contrasts are also sensitive to as shown in the panel on the right. 



D. Example III: Q = 




Dark-sector interaction of the form 



Q 



(70) 



has been previously analysed in [3, ll3|, l33|, |36|, |37|, |33| . This type of interaction appears in some scalar-tensor theories 
with an interaction vector of the form 



(71) 



where T^"^) = —pc is the trace of Tj^l! . Here, we loosely refer to (f> as quintessence although it is strictly more than just 
another form of cosmic energy density. In a scalar-tensor theory, (f> is a mediator of gravitational interactions and can 
only be interpreted as part of the dark sector after a conformal transformation to the Einstein frame. 

By perturbing the vector ([TTjl and following the procedure in the previous section, we find the perturbed interaction 
to be 



6Q 



ic) _ 



(72) 



In terms of the dimensionless variables, the interaction and its perturbation are 

7 — bu^ , 



5^ = hu^ Sc + -^h{K.5(j}) 

V v6 

Substituting the above into ((24| -(|25 p gives the coupled differential equations governing the perturbations 



(73) 
(74) 



dN^ 



dlnH , \ dSc 

\-2-2bx] — 

dN / dN 



iVebx^ i 1 



26^ - ^ ] u^6c = 2\/6 



46 



b b^x 

- + x 

3 3 



d{K5(j)) 
dN 



2V6b 



dlnH 16 



dN V6 



^b^x - ^(-^] + 2V6b - Vebrjy^ + 3Xy'\2b^ + 1 
v6 v6 \aH 



dlnH 



dN 
dSc 



dN 



3 + 4bx 
~ Vdbu^Sc 



d{K54>) 
dN 



3^y^ + ( — ] + 6&^i^ + 46^2;^ - 2V6bXy- 
' all ' 



(kScI)), 



b^u' 



(75) 



(76) 



Note that there is no longer any singularity near x = 0. 



1. At the scaling solution 



For this model, the scaling solution exists whenever 



'-^<Bb<^{B--,)^ (77) 
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FIG. 4: Trajectories in the x — y plane for the exponential potential with A — \/5 and interaction term Q — ^j2j?ihK(\)pc. From 
left to right, the trajectories within the semicircle correspond to 6 = 1, 0.5, 0, —0.5 and —1. The dotted hyperbola is the locus 
of scaling solutions (Equation (I61|l 'l. In contrast with Figure ((T|, the region a; < is now accessible to the trajectories since 
there is no singularity at a: = 0. 




FIG. 5: Evolution of large-scale pe rtur bations over AA^ = 10 during the scaling regime for the exponential potential with 
A = %/5 and interaction term Q — ^2/3tiK0pc (with & = 0, ±0.1). The panel on the left shows the evolution of density contrast, 
5c, while the panel on the right shows the quintessence energy contrast, The dashed lines indicate negative values. 



where B = ^5 + A [H, with 

_ 3 _ \J^\? + 2\/6A& + 9 _ v/6A2 + 2V6A6 - 18 dlnF _ Sv^A 

^ ~ %/6A + 26' ^~ V6A + 26 ' "~ V6A + 26 ' dN 2(\/6A + 26)' ^ 

Surprisingly, the locus of the scaling solutions remains the same hyperbola (|6ip . The background dynamics in the 
phase space (x, y) is shown in Figure [4] Note that trajectories are now free to explore the entire semi-circle. If 6 > 0, 
this type of interaction gives rise to an effective opposing force that slows down the field, just as in the previous case. 
In this case, the field can even roll 'uphill' without encountering any instability. Similarly, we see that interaction 
with 6 < boosts the velocity of the field. 

Figure [5] shows the results of numerical integrations of the density contrasts be and 8^. These results are similar 
to the previous case, i.e. the growth rates of the perturbations increase with 6 < while the reverse occurs with 
6 > 0. The gradients of the lines representing decaying and growing mode are also approximately —3/2 and 1 as 
in the previous case. However, the growth rates in this case respond more sensitively to changes in the interaction 
parameter. 
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FIG. 6: Evolution of large-scale perturbations over AA^ = 10 during the quintessence-dominated regime for a flat potential 
V ^constant (A = 0) and interaction term Q = y^2/3bK,<j>pc. The panel on the left shows the evolution of density contrast \5c\. 
Values of b up to ±0.5 are shown for comparison. The panel on the right shows the quintessence density contrast \S^\. The 
asymptotic forms of these curves are discussed in the text. 



2. Flat limit 



The position and stability conditions for the quintessence-dominated attractor remain the same as in the previous 
coupling. Substituting in the coordinates (|49p and using the large-scale and flat limits, the perturbation equations 
([75)) - ((76)) can be solved exactly, with solutions 

5c = Ci + C2N + Cse-^'^ , KS(t> = + de-^^ , (79) 

where Ci are constants and 6 7^ 0. The case b = was covered in section HVBI Interestingly, there is a new linear 
mode in Sc, and, unlike the previous case, the quintessence-induced mode in Scj) is simply constant. This implies that 
Sc evolve linearly with N, while remains constant. This is in contrast with the case without interaction, in which 
decays exponentially to zero. 

Figure [S] illustrates the features outlined above. In the panel on the left, in which a linear scale is used, the adiabatic 
condition induces a linearly decaying 6c for 6^0, whilst 6c remains constant for 6 = 0. The difference between the 
two cases are only noticeable with larger values of the interaction, shown here up to 6 = ±0.5. 

The panel on the right of Figure [S] shows the evolution of the large-scale quintessence perturbation 6^. The constant 
modes are clearly seen in the cases where b ^ 0. Without any interaction, quintessence perturbations decay with slope 
—3, as expected from Equation ((50|) . Over AiV = 0, the difference between the coupled and uncoupled cases grows 
to some 20 orders of magnitude. Consequently, the details of structure formation and large-scale CMB anisotropics 
in these cases will also be markedly different. 

Moreover, such a vast difference in the magnitudes of dark-energy perturbations persists down to small scales, as 
shown in Figure [7] (in which k > aH in —10 < iV < 0). Without interaction, (5^ decays to zero as before. For 
b = ±1, 1 5^1 is no longer constant on small scales, but oscillates briefly before a growing mode dominates^. This 
makes it possible for dark energy to clump on sufficiently small scales to have important astrophysical consequences, 
for instance, on the abundance of galaxy clusters and on gravitational lensing. 



V. CONCLUSIONS 



In this paper, we have investigated the consequences of interactions between dark matter and dark energy in the 
form of quintessence. We have set up a general framework in which the density contrasts in dark matter and dark 
energy could be calculated given any quintessence potential V{(j)) and dark-sector interaction Q^. Our formalism 



^ We find that Sc is also enhanced on small scales, in agreement with the results of Il9ll . However, such large perturbations should be 
confirmed by a calculation which uses a nonlinear perturbation theory. 
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N 

FIG. 7: Evolution of small-scale quintessence density contrast during the quintessence-dominated regime for a flat potential 
V ^constant (A — 0) and interaction term Q = ^2/3&K(^pc. The dashed portions indicate negative values of 5^. 



is built on the phase-space method in which different cosmological histories correspond to trajectories in a finite 
space of dimensionless variables (Pt- pH]) . We have shown how the phase-space formalism may be extended so that 
perturbations in the densities of dark matter and dark energy may be calculated along the trajectories. This is 
particularly useful when combined with the knowledge of the attractor dynamics in the phase space. 
The key results of this paper are as follows: 

• Using the results of Hwang and Noh [llj , we obtained a set of second-order coupled differential equations ([M]) - 
(j25p describing the evolution of the perturbations in terms of dimensionless background variables [Equation 
(j26p ]. This is valid for any quintessence potential and interaction vector. We note that the perturbation in 
the interaction itself plays a crucial role as extra source terms in the coupled system. 

• We explained why adiabaticity is no longer conserved in the presence of interactions, even for large-scale per- 
turbations [Equation (^1)) ]. 

• Our formalism was applied to a number of scenarios involving the exponential potential. In particular, we 
studied the behaviour of perturbations as trajectories approach the scaling solution as well as the quintessence- 
dominated attractor in the case that the potential is almost flat. To test this technique, we set the interaction, 
Q, to zero and obtained analytic results that agree with those obtained by previous authors. 

• We considered the interaction of the form Q — (3Hpc. At the background level, we found that /3 > effectively 
acts against the motion of the field, reducing its kinetic energy. Applying the perturbation equations, we found 
that the growth rates of perturbations change very slightly in the scaling regime. At the quintessence-dominated 
regime, we found analytic expressions for the growing and decaying modes for large-scale perturbations. 

• An instability associated with the interaction Q = PHpc was found. When the field velocity is close to zero, the 
perturbation equations blow up despite the fact that the background variables continue to be well-behaved. We 
believe that this is the first time that this type of instability has been identified. To avoid such an instability, 
the initial conditions must be tuned so that the field consistently rolls without stopping. To the author, such 
an extra tuning of initial conditions spoils the dynamical advantages of quintessence and makes this type of 
interaction unappealing. 

• We investigated the coupling Q = •\/2/3&K^pc- We found that b > works against the motion of the field just as 
before. However, perturbative instabilities are absent in this model even when the field rolls 'uphill'. In general, 
we conjecture that all models of dark-sector interaction with Q — 0{(j)^^) and n < 1, suffer from instability in 
the perturbations when the quintessence velocity is small. 

• For 7^ 0, dark-energy perturbations remain constant in the quintessence-dominated era instead of decaying 
away, even on small scales. Figure [7] shows an example in which the difference in the dark-energy perturbations 
between the cases 6 = and |6| — 0.1 grows to over 20 orders of magnitude. This opens the possibility that dark 
energy could clump on astrophysical scales. The implications for the CMB anisotropies, matter power spectra, 
galaxy cluster counts and gravitational Icnsing observations are interesting, and we shall investigate further into 
these issues. 
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As long as there is no conclusive evidence otherwise, one must remain open to the possibility that there may be 
interactions between the two most enigmatic components of the Universe. Nevertheless, as we have shown in this 
work, the mathematical details of dark-sector interactions can be challenging and must be handled with care. There 
are good prospects of identifying such interactions since a small interaction coefficient may give rise to a dramatic 
increase in the magnitudes of dark-energy perturbations. It would be interesting and quite straightforward to apply 
our formalism to other types of potentials and interactions. 
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